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Abstract 

A discrete A^— level alternative to the popular imaginary cubic oscillator is 
proposed and studied. As usual, the unitarity of evolution is guaranteed by 
the introduction of an ad hoc, Hamiltonian-dependent inner-product metric, 
the use of which defines the physical Hilbert space and renders the Hamil- 
tonian (with real spectrum) observable. Due to the simplicity of our model 
of dynamics the construction of the set of eligible metrics is shown tractable 
by non-numerical means which combine the computer-assisted algebra with 
the extrapolation and/or perturbation techniques. 
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1 Introduction 



Many measured spectra of energies may be interpreted as excitations of a 
quasiparticle. The simplest fits of such a type employ the elementary one- 
dimensional differential Schrodinger equation 

(P 

-^^n{x) + V{x)ij{x) = EniJn{x) ^(±A) = , A < oo (1) 

containing a real potential V{x). In more sophisticated models, potential 
V{x) may even be allowed complex, provided only that the spectrum itself 
remains real (cf. reviews [H El |3l H] for details). 

Whenever V{x) ^ V*{x)^ the reality of the spectrum may be fragile 
and sensitive to perturbations [5]. A remarkable exception emerges with 
the robustly real spectra generated by many potentials with the property 
V{x) = V*{—x) called, in the hterature, PT— symmetry ^ alias parity- 
pseudo-Hermiticity [3] alias Krein-space-Hermiticity [B]. 

In the early studies of this remarkable mathematical phenomenon the 
attention of the authors has mainly been restricted to the imaginary cubic 
oscillator example V{x) = V^^^\x) = ix^ modified, possibly, by some other, 
asymptotically subdominant terms (cf., e.g., papers by Caliceti et al [7], by 
Alvarez [|8J or by Bender et al [9l and several further authors [H]). In 
what follows we shall mainly feel inspired by this choice of V{x) as well. 

The essence of our present message will lie in the recommendation of a 
drastic simplification of the necessary mathematics. In a way explained in 
section [2] this will be achieved by means of the replacement of the differential 
Schrodinger equation ([1]) by its discrete, difference-equation analogue. A few 
simplest illustrations will be then added in section [3] where the discretization 
of the coordinate will be shown to facilitate the study of the reality (i.e., in 
principle, observability) of the spectrum. 

The desirable flexibility of our present discrete simulation of the dynami- 
cal energy spectra will be shown achieved by a supplementary one-parametric 
deformation of the potential resembling slightly the influential proposal of de- 
formation ix^ — )• {ix)^~^^ by Bender and Boettcher |10j. The details will be 
described in section |H The subsequent section |5] will then recall the known 
theory and explain some of its details via the simplest possible example with 
N = 2. In sections |6]-|8] we shall finally present some applications of this the- 
ory to the models with = 4, = 6 and general A^ = 2K > 8, respectively. 
Section M is summary. 

2 Discrete Schrodinger equations 

The spectrum of many Krein-space-Hermitian Hamiltonians H ^ W has 
been found robustly real and bounded below [101 [121 [IS]- In the spirit of the 
general theory as outlined, first, by Scholtz et al [H], one can conclude that 
the apparent non-Hermiticity is "false" and that it may be reinterpreted as 
a mere consequence of an inappropriate choice of the Hilbert-space repre- 
sentation L^(]R) = 'H}~^\ Hence, the abstract remedy is straightforward and 



lies in an interpretation-mediating transition to a "standard" Hilbert space 
of states 7^(^) [3]. 

In the practical applications of such a theoretical scheme the original 
vector space is usually being endowed with a general, non-/oca/ inner product, 



= / / ^*{x)Q{x,y)4>{y)dxdy={ij\e\^)^^\ 6 = 0^ > . (2) 




For the one-dimensional imaginary cubic oscillator Hamiltonians 

(where the dots may represent certain asymptotically subdominant terms) 
the first constructions of the necessary "standard metric" were perturbative. 
They appeared in 2003 [13] and in 2004 [16j, yielding the metric operator 

e = e(//(^^)) = exp(-giA - QsA' - . . .) (3) 

with Qi{x,y) = ixy{x^ + y^)sign{x — y)/2'^, etc. In the light of the well 
known fact that the assignment of the metric G to a given Hamiltonian H 
cannot be unique in general [T^ [T7] . additional, multiparametric metrics 
0(ci,c2,...)(-^''^'"^) were further found and constructed in Refs. [16] and [18]. 

The latter calculations proved facilitated by an additional assumption of 
the absence of a nontrivial fundamental length L > in the theory. Un- 
fortunately, the presence of such a preassigned length scale has been found 
essential in Ref. [19] where an extension of the formalism to the scattering 
dynamical regime has been proposed. Naturally, under the assumption of 
the presence of a fixed length scale L > 0, additional ambiguities will emerge 
in the metrics. In the differential-operator models their explicit specification 
might prove prohibitively difficult (cf., e.g., Ref. |i20j or section 6 of Ref. [IB] 
for related comments). 

We shall address here this problem while accepting the discretization 
strategy of Ref. [I9]. In the way based on the use of equidistant, Runge- 
Kutta grid-point coordinates Xk = —A+k h , k = 0,1, . . . , A^+1, the ordinary 
differential Schroedinger Eq. ([T]) will be replaced by its discrete version 

_ ^fa-0-2^fa) + ,i.fa,.) ^^^^^^ ^^^^^ ^ ^ ^^^^j 

with xn+i = A (i.e., h = 2A/{N + 1)) and with the Dirichlet boundary 
conditions 

ipixo) = ijixN+i) = . (5) 
This leads to the reduction of the problem to the (numerical) determination 
of the A^— plets of eigenvalues Ej^\ j = 1,2, . . . , N , i.e., of the eigenvalues 
Ej := h'^Ej^'^ G (0,4) of the re-scaled, A^— dimensional Hamiltonians 

f2 + h^Vixi) -1 \ 

-1 2 + h^V{x2) -1 

i/(^)= -1 2 + h^V{xs) ■•. 



2 + h'V{xN) I 



(6) 



We intend to make use of this definition of the Hamiltonian in the whole 
rest of our present paper. Nevertheless, before we fully concentrate on its 
bound-state aspects and consequences, let us make a small detour and point 
out that our present acceptance of the Dirichlet boundary conditions is 
in fact just one of the two basic alternative options which are at our disposal 
in the phenomenologically oriented applications. In a very close parallel to 
the continuous cases, these boundary conditions could have been replaced, 
even in the Runge-Kutta-discretized models, by the alternative scattering 
scenario. 

In the latter context, one might feel discouraged by a threatening tech- 
nical difficulty connected with the necessity of working with the infinite- 
dimensional matrices. Still, there exist tricks which enable us to avoid these 
difficulties by assuming that our interaction V{x) in Eq. is just short- 
ranged (plus, admissibly, non-Hermitian and non-local). For more details 
and/or for an explicit constructive illustration of such an alternative possi- 
bility, interested readers might consult, e.g., Ref. fIT\ . 



3 The series of one-parametric Hamiltonians 



Let us pick up the potential V{xj) = ix^ and insert it in Eq. 



Once 



we abbreviate a = h^/8 we obtain the following one-parametric sequence of 
Hamiltonians 
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-1 2 + ia 
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(7) 






-1 
2 + 27ia / 

of this series is particularly elementary. For this 
reason it has also been chosen and studied in the methodical part of Ref. [15] 
(cf. its section II. B). 

For the sake of brevity we shall restrict our attention to the matrices 
with even dimension N = 2K. Thus, in the simplest case with K = 1 
one encounters the compact and explicit energy formula £1,2 = 2 =]= yl^^c?. 
Such a two- level spectrum remains real iff a G (—1,1), forming a circle in 
the energy-coupling plane. The similar formulae and conclusions may be also 
obtained at the next few Ks. 

The first nontrivial sample of the a— dependence of the spectrum may 
be obtained at = 4. Its shape is displayed in Fig. [TJ The algebraic 
representation of these energies is elementary. 



2,3,4 = 2 + 1/2 V 6 - 1460 a2 ± 2 V529984 a^ - 1680 a^ + 5 . 



One can easily deduce that these energies remain real for a G (- 
where a^^) = 1/2 - 1/18 ^69 ^ 0.0385208965. 
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Figure 1: The parameter-dependence of the real-energy eigenvalues e{a) of 
the toy Hamiltonian H^^\a) [= the third item in the hst ([7])]. In topological 
language, this spectral locus is formed of the two (deformed) cocentric circles. 
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Figure 2: Same as Figure [H, at the next even dimension = 6. The spec- 
tral locus is formed of the three (deformed) cocentric circles (note that the 
horizontal axis is rescaled). 

Numerical difficulties start emerging at = 6 since the purely alge- 
braic representation of the spectrum requires the use of Cardano formulae in 
which the imaginary numerical errors occur and survive, disappearing only 
in the infinite-precision arithmetics. Still, the more or less routine control of 
precision enables us to conclude that the N = 6 spectrum remains real iff 
a e (^—a^^\ a^^^Y The approximate numerical value of a^^^ ~ 0.011344897 
may be read out of Fig. H] and/or of its appropriate systematic magnifications. 

We may notice that at = 4 and = 6 the ends ia*^^'^-' of the interval of 
the reality of the whole spectrum are determined by the confiuence (followed 
by the complexification) of a single pair of energies (lying in the middle of the 
spectrum). Both of these models are exceptional. Starting from N = 8 the 
pattern gets changed and we encounter the less trivial a— dependence of the 
energies e{a) characterized by the fragility and confiuence of the two pairs of 
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Figure 3: The next item in the series of Figures [T] and O The horizontal 
axis is further rescaled and the = 8 spectral locus gets nontrivial. In 
topological terms, it becomes composed of the vertical array of the three 
circles, all being circumscribed by the fourth one. 

Table 1: The exceptional-point coordinates of the loss of reality of the first 
and second excited state {z = 3). 



dimension 


critical 


twice-degenerate 




parameter a 


energy £2 (a) = £3(0^) 


2 


1 


2 


4 


1/2 - V69/18 


2 


6 


0.011344897 


2 


8 


0.003383828 


0.9285 


10 


0.001246890 


0.6194 


12 


0.000543788 


0.4438 


14 


0.000266880 


0.3335 



the energy levels (cf. Figj3]). 

The latter form of the loss of reality is generic. The = 14 sample of 
the spectrum as displayed in Fig. H] elucidates some details. In the picture we 
see that the spectral locus (i.e., the set of all of the real energies £j{a) with 
a G (— cxD, cxo) and j = 1,2,..., N) preserves the form which remains up- 
down and left-right symmetric. During the growth of = 2K the graphical 
analysis of the model remains feasible and reveals a steady shrinking of the 
interval of the parameters a for which the spectrum remains real. 

A quantitative account of the latter phenomenon is presented in Table [H 
The inspection of this Table reveals that the presence of the fundamental- 
length parameter a^^^ (= a maximum of admissible as, decreasing with A^) 
in our present model implies an obvious mismatch between the large— A^ 
spectrum (with a very small interval of reality) and the robustly real spectrum 
which has been proved to exist in the differential-equation limit N = 00 [12]. 

An explanation of the paradox is twofold. Firstly, the enormously small 




Figure 4: A sample of a "multilevel descendant" of the preceding Figures. 

At = 14 we see a "generic" pattern in which we encounter (A^ — 6)/2 
concentric circles plus a single "top" and single "bottom" circle, all being 
circumscribed by the "biggest" last circle. The interval of the reality of all 
of the eigenvalues e{a) of the toy Hamiltonian H^^'^\a) has further shrunk. 

magnitude of a — 0{h^) makes the comparison only sensible in the limit 
a ^ (in this sense, the spectrum of the discrete model remains robustly 
real at all A^). Secondly, at any finite dimension N < oo even the role of 
the Runge-Kutta error terms ^ 0{a) remains unspecified. Whenever 

the limit of A^ — )■ cxo is considered, this observation offers another argument 
against drawing any A" — > oo implications from the observations made at the 
finite A" and nonvanishing values of a in our models. 

4 The series of two-parametric Hamiltonians 

In the examples of the preceding section we encounter another suspicious 
feature even at any fixed and finite dimension A^ = 2K. Indeed, the critical 
left and right points a — ±a^^^^ of the loss of the reahty of the spectrum 
seem to be exclusively connected with the confluence of the flrst and second 
lowest excited states £2 and £3 or, symmetrically, of their equally "privileged" 
mirror partners eN-2 and En-i- 

In what follows we shall explain this apparent privilege as an artifact 
caused by the too speciflc choice of the potential. The constructive expla- 
nation will be based on a rescahng of the interaction V{x) — ix^. Its third 
power term will be replaced by an odd function f{z){x) ~ sign(a;) with a 
real exponent 2; G M. As long as our models are discrete, we need not follow 
the conventional wisdom and require that the functions f{z){x) are analytic. 



•7 



Thus, we arrive at the family of two-parametric Hamiltonians H^'^^\a, z) = 
-2-ia{2K-iy -1 ... ... 

■-. 2-za3^ -1 

■ ■■. -1 2-m -1 

-1 2 + m -1 ■•. : 

-1 2 + ia3' '■■ 

... ... -1 2 + ia{2K-iy 

(8) 

We shall see below that such a class of discrete deformations 

V{x) =ix^ := if(3){x) V{x)= if(z){x) 

may be made responsible for the emergence of a full range of possible com- 
plexification patterns. In the language of dynamics and phenomenology, we 
shall show below that the emergence of the new parameter z in models ([8]) 
renders them tractable as offering a certain complete set of alternative, tune- 
able topological features of the spectrum. In this sense our present complex- 
interaction models may be perceived as complementing the recent attempts 
of achieving a topology-related tuning of spectra via a non-locality of real 
V{x) [22] or, alternatively, via a nontriviality of the "coordinates" x living, 
say, on real closed loops [23] or on certain specific multisheeted complex 
curves 

4.1 Spectral loci e^^'^)(a): the = 4 prototype 

As we indicated, the key source of interest in our two-parametric mod- 
els H^'^^\a, z) should be seen in the enhancement of the flexibility of the 
a— dependence of the spectrum. Moreover, one can expect that the choice of 
z might re-assign the role of the most fragile levels all along the spectrum. 

In an expected verification of these hypotheses, let us now select = 4 
and turn attention to the two-parametric Hamiltonian 

2-ia3' -1 
-1 2-ia -1 

(9) 

-1 2 + m -1 

0-1 2 + ia3^ _ 

Due to the simplicity of such a generalization of the purely cubic discrete 
model ([7]) it may be easily shown, graphically, that the spectral pattern 
displayed in Fig. [T] does not change too much when the exponent-parameter 
z starts to be different from three. Indeed, with the steady growth of z > 3 
one merely reveals that the small internal ellipse in the picture of Fig. [1] 



H^^\a,z) = 



will shrink. In a numerical test performed at z = 25, for example, we still 
recognized the existence of this inner ellipse but only on the scale of as of 
the order of magnitude of 10~^^. 
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Figure 5: The a— dependence of the real eigenvalues e{a) = e^^'^\a) of the 
two-parametric toy Hamiltonian H^^\a,z) ai N = 4 [cf. Eq. (^]. Slightly 
above a critical value of z = 85/64 ~ 1.328125 > ZcnUcai, the pattern is still 
formed by the two cocentric circles. 

In the opposite direction of the change, i.e., with the decrease of z below 
three the internal ellipse broadens. From the inspection of Fig. [5] we may 
deduce that slightly below the value of z = 85/64 the internal ellipse must 
ultimately touch the external deformed circle. At this moment the whole 
spectral locus will acquire the form of the two intersecting ellipses. After the 
further small decrease of z we recognize a qualitatively (i.e., topologically) 
new pattern, the form of which is sampled, at the next rational numerical 
value of z = 84/64, in Fig. [HI 
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Figure 6: The breakdown of the cocentric-circles pattern of Figure [5] when 
the value of z dropped just slightly below the critical value, z = 84/64 = 
1.3125 ^ Zf^j-iiiQui. 

With the continuing decrease of z we witness a quick shrinking of the 
two spurious large- |a| intervals of partial reality of the spectrum. Both of 



n 



them are centered around \a\ = 1/2 while leaving still the two levels real 
(cf. Fig. [7] where we choose z = 81/64 = 1.265625). We also repeated the 
same graphical analysis below the latter value of the exponent. We revealed 
that quickly, both the anomalous partial-reality intervals have got empty and 
disappeared. 
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Figure 7: The "evolution-to-disappearance" of the left and right "anomalous" 
topological circles of Figure |6] at 2; = 81/64 = 1.265625 < ZcHUcai- During the 
further decrease of 2;, just the real spectral locus composed of the horizontal 
array of two nonintersecting (deformed) circles will survive. 

Needless to add that the further steady decrease of z already keeps the 
"two-oval" topology of the spectral pattern unchanged. We checked this 
empirical rule in the exactly tractable limit of 2; = 0. In this limit the 
quadruplet of the energies acquires the compact form 




One can strictly deduce that this spectrum remains real iff a G (-^5/4,75/4). 



4.2 Sequences of rearrangements: the = 8 prototype 

The straightforward generalization of Eq. ([9]) to any dimension yields the 
N hy N matrix H^'^^\a, z) of Eq. (j8]). The related geometric shapes of the 
spectra just generalize the = 4 pattern. At any N = 2K and at all 
of the sufficiently large exponents we revealed that the picture of the real 
spectral locus in the (e, a)— plane remains topologically equivalent to the set 
of cocentric circles. At = 8 and 2 = 9/2 and z = 8/2, the respective 
Figs. [8] and |9] offer the two characteristic samples of the evolution of the 
large-exponent spectrum with the decrease of the exponent z > z^j^^J^ critical- 
In a series of graphical experiments using the smaller and smaller rational 
exponents z we were able to keep the numerical precision under good control. 
We discovered that Zj--l^^ critical ~ 4- Below this value though safely above 

^second critical ~ ^2/4 (cf. Fig. [3]above), one encounters the new topological 
pattern sampled by Figs. [10] or [TTl 
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Figure 8: For the less elementary matrix H^'^^\a, z) of Eq. ([8]) with K = 
the spectral locus remains composed of the four concentric (deformed) circles 
at all the sufficiently large exponents z. This pattern is sampled here at 

^ ~ 9/2 > ^ fir St critical- 
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Figure 9: The topology of Figure [HI is not yet changed in the "close-to- 
crossing" case oiz = 8/2 > zfl^^ 



The corresponding second-pattern regime may further be split in two 
stages. In the initial stage of the decrease of z there survive two large— |a| 
anomalies which quickly shrink and, not too far below the first critical value 
of the exponent, disappear. In the final stage of the decrease the picture of 
the whole real spectral locus in the (e, a)— plane acquires the "fully canonical" 
form of the two cocentric (deformed) circles with a vertically ordered pair of 
non-intersecting circles (or rather deformed ellipses) inside. 

The next, second change of the topological pattern has been spotted to 
occur between z = 13/4 (cf. Fig. [TTi) and z = 7/4 (cf. Fig. [T^D while the third 
change certainly follows between z = 7/4 and z = 6/4 (cf. Fig. [T^ . Ulti- 
mately, the list of the topological changes of the N = 8 spectral loci is made 
complete during the transition between z = 6/4 and z = 1/2 (cf. Fig. UM . 
One arrives at the other, small-exponent extreme in which the spectral locus 
stays unchanged and topologically equivalent to a vertically ordered quadru- 
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Figure 10: The innermost circles of Figures [H] or |9] are being replaced by 
their horizontal array ai z = 14/4 G (^^^^ critical^ ^Lt critical)- Notice that 

(8) 

we moved safely below Zj/j.^-). critical ^o that the temporary left and right real- 
eigenvalue anomalies have already (and very quickly) disappeared. 

plet of separate deformed circles. 

The results of our calculations aimed at a more precise determination of 
the numerical values of the quadruplet of critical exponents at = 8 are 
summarized in Table O In order to make these calculations feasible, we were 
forced to give up the full control of the numerical precision. For this reason, 
Table [2] does not contain any estimates of the error bars so that even the last 
digits in our numerical values of the critical exponents z should not be taken 
for granted. 

Marginally, let us add that, strictly speaking, our present list of the five 
critical points might have been complemented not only by the list of cer- 
tain "secondary critical" points z (marking the disappearance of the above- 
mentioned left and right partial-reality intervals in the anomalous, large— |a| 
dynamical regime) but also by the list of certain "tertiary critical" points z 
at which these left and right anomalies change from "containing" to "not- 
containing" an even narrower four-level-reality subinterval. 

The readers who would be interested in the similar subtleties might con- 
sult the results of this type as obtained, e.g., in Ref. [22]. In the present 
context, a sample of such an analysis has merely been performed in the 
small-exponent regime with z > where we found z^^^^rth critical ~ 1-0358. 

In this setting the unchanged topology without anomalies has been only 
demonstrated to exist in the slightly smaller interval of z & (~oo, z) with 
z ~ 1.033. Within the latter, anomaly-free interval, the variations of z still 
preserved the strict reality just along the four non-intersecting ellipses, more 
and more deformed and all the time located in a vertical arrangement. In 
the adjacent, very short interval of 2; G (1.033, 1.0358) there emerged the left 
and right anomaly as sampled, at = 4, in Fig. [71 

At the fourth critical value of z ~ 1.0358 we obtained, in accord with our 
expectations, the transitional intersection of ellipses as sampled, at = 4, 
in Fig. [6] We can conclude that our numerical study of the = 8 model 
confirmed that when we ignore the partial-reality anomalies as inessential 
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Figure 11: The topology-preserving deformation of Figure [TO] during the 
further decrease of the exponent to 2; = 13/4 > z^J^^^^ critical (notice also the 
change of scale of a). 




Figure 12: The next change of topology of the spectral locus as sampled at 

^ ~ '^/^ ^ ^third critical- 



for the fully unitary quantum systems, the range of the exponents z > 
(or rather z > —00) splits into five subintervals on which the spectral loci 
e{a) become topologically non-equivalent. The exhaustive and fully explicit 
description of these N = 8 topologies is also given in Table El 

4.3 Rearrangements: Fibonacci-sequence connection 

Fig. [8] may be perceived as the inessential N = 8 sophistication of the = 4 
pattern of Fig. [5l Similarly, Fig. [9] recycles the structure which is shown in 
Fig. ini while Fig. [10] may be read as paralleling Fig. [71 The generic > 6 
pattern only starts emerging, at = 8, during the transition from Fig. [Til 
to Fig. [T2I It is characterized by the two parallel confluences of energy pairs 
at the critical exponent Zsecond critical- The similar phenomenon occurs at all 
of the further critical points, viz., at the remaining two topology-changing 
transitions from Fig. [121 to Fig. [131 and from Fig. [131 to Fig. [TH 
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Figure 13: The next topological arrangement of the spectral locus sampled 

at ^ = 6/4 G {z^f^-^^ critical y ^fourth critical)- 
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Figure 14: The ultimate, "small-exponent" topology which is sampled here 
at 2 = 1/2 but which is exhibited by any spectral locus e{a) of matrix ([8]) at 
= 4 and any real z < zf]^^^^ 



At any higher dimension = 27^ the similar sequence of changes of 
topology can be assigned a fully systematic description. Thus, once we choose 
any and assume the complete knowledge of the set T/v of all of the related 
spectral-locus arrangements (consisting of the /T— plets of deformed circles) 
which are topologically non-equivalent, we may contemplate a constructive 
transition to the next dimension A^' = IK + 2. In the first step we build the 
first subset T^i of new schemes T/v' as the set of elements of the old set T/v 
which are merely encircled by an additional, single outer (deformed) circle. 

The construction of T/v' = ^n' S^^i may be then completed by the 
discovery of the one-to-one correspondence of the second part S^i of the new 
set of schemes to the "one step older" set Tat" where A^" = 2K — 2. Each of 
the "older" elements only has to be complemented by the jtair of the single 
upper and the single lower additional (deformed) circles, indeed. 

In the enumeration of the complete sets of non-equivalent patterns by 
induction, it is now sufficient to verify that the numbers j^T2K of the elements 
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Table 2: The five intervals of exponents yielding topologically non-equivalent 
real spectral loci £:(a) at = 8. 



the graph of ej = Sj{a), j = 1,2, ... ,8 


illustration 


interval 


four cocentric (deformed) circles 
two in vertical array, encircled by two 
three in vertical array, encircled by one 
vertical, but the middle pair cocentric 
four (deformed) circles in vertical array 


Figs. El i 

Figs, m M 

Fig.m 
Fig. [13] 
Fig.m 


z e (3.982, oo) 
z e (3.178,3.982) 
z e (1.630,3.178) 
z e (1.0358,1.630) 
z G (-CX), 1.0358) 



of the sets T2K are equal to K for K = 1, 2, 3. Thus, we have proved 

Lemma 1. We have #72ic = F^-i where Fj is the j — th Fibonacci number. 

Remark 1. The sequence of Fibonacci numbers Fj is defined by recurrences 
Fj = + F,_2 yielding i^Ts = 5, #rio = 8, #ri2 = 13, etc. 

5 Interpretation 

Usually [25] people decide to work in a fixed, specific Hilbert space T-L^^^ and 
treat a given quantum system as physical if and only if its time evolution is 
unitary. Naturally, this is the strategy which does not change if the space 
proves endowed with a general, nontrivial metric 6 = 0^ > (i.e., with 
its inner product defined in terms of this metric, cf. Eq. ([2]) above). One still 
has to guarantee that every candidate for an operator of observable proves 
self-adjoint with respect to this metric [14\. 

For our present, manifestly non-Hermitian matrix representations H^'^^ 
of the Hamiltonians (which will have to play the role of the generators of the 
unitary time evolution in 'H^^'') we must guarantee their Hermiticity with 
respect to the given nontrivial metric. Such a form of nontrivial Hermiticity 
could be better called cryptohermiticity In the context of mathematics, 
such a condition is often being interpreted as the Dieudonne's [26] "quasi- 
Hermiticity" constraint. 

In the application of such an approach which is to be employed in what 
follows, we shall always start from a given matrix H^'^^ and reconstruct the 
ad hoc metric (or metrics) 9 = 9 (if) via Eq. ([TOl) . treating this constraint 
as a linear algebraic set of equations for the matrix elements of the metric. 

5.1 The simplest N = 2 illustration 

For the most elementary, z— independent N = 2 "input" Hamiltonian (cf. 
the first item in Eq. ([7])), the solution of Dieudonne's constraint f[TO|) yields 
the complete, two-parametric family of the most general matrices of "output" 



1 c 



pseudometrics, 



(k,m) 



k km — ika 
km + ika k 



k,m E 



(11) 

Physical condition a G (—1,1) guaranteeing the reahty of the energies must 
be complemented by the condition of positivity of the (two) eigenvalues 



= 0^ = k± -Jk^m^ + k^a^ > . 



(12) 



Only such a condition will open the possibility of using the corresponding 
pseudometric ( fTTj) metric in Ti^^^ [T4] . 

We see that k must be positive and larger than the square root. We may 
set a = cos/3 sin 7 and m = cos/3 cos 7 with, say, /3 G (0, vr) and 7 G (0,7r) 
and — 1 < cos /3 < 1 . This reparametrization leads to the final and most 
general positive definite and Hermitian metric 



e = e{i:f(2)[a(/3,7)]} 



[fe,m(/3,7)] 



k- 



1 

'-'P COS /3 



-17 



cos /3 



(13) 



which is never diagonal. As long as it contains free parameters, their choice 
will fix the inner product and, hence, it will specify the Hilbert space 
of admissible states of the model. The choice of parameters /3 and 7 will 
determine both the Hamiltonian H^'^\a) and the metric 0|f^„) (with, say, 
/c = 1), i.e., the true and complete physical contents of the theory. 



5.2 Eligible observables 



(2) 

In our N = 2 theory, any matrix A = A^^\^-| representing an observable 

lis) 



quantity must be self-adjoint in w^iJ^y i.e. we must have 



A{^^^) 0[fc,m{/3,7)] - 



[fc,m(/3,7)] A(/3,^) . 



This equation must be satisfied by the eligible representation matrices 

G + ig B + ib 



A 



C + ic D + id 



(14) 



(15) 



Naturally, Eq. ffT^ admits an arbitrary A;— rescaling of G |if*^^^[a(/3,7)]} 
and/or a similar trivial rescaling of A(/3 .y). In terms of matrix elements one 
can easily check that this equation imposes four real constraints upon the 
eight free parameters in ( IT5l) . This means that the N = 2 family of available 
observables is four-parametric in general. 

The details of the construction are left to the readers. We can only 
summarize that the three constraints are trivial and that they merely define 
quantities B, C and the difference G — D. The remaining, fourth constraint 
acquires the form of a linear relation between sums Ce = 6 + c and Qt, = g + d 
with the unique solution gj] = 0. One can conclude that for a given input 



m = m(/3, 7) and a = a(/5, 7) the final and entirely general form of the N = 2 
observable reads 



A = AiD,b,c,g) = - 
a 



Da — b — c + iga, g — bm + iba 
g + cm + ica, Da — iga 



(16) 



I.e., 



We can check that our original Hamiltonian is reobtained at D = 
G = 2),b = c = (i.e., B = C = -1) and g = -a (= -d). 

In Ref. [13] the authors recommended to proceed in an opposite direction 
and select a few matrices (or operators) Ai, A2, . . . of observables in advance 
(say, on some fitting or phenomenological grounds). Naturally, these input 
matrices must necessarily be self-adjoint in 7i^^\ i.e., in our N = 2 example, 
proportional to our general formula f lTB]) . Thus, in the generic case, the series 
of the necessary cryptohermiticity conditions will, sooner or later, specify all 
of the values of the free parameters in the metric (up to the above-mentioned 
trivial rescaling of course). 



5.3 The concept of charge 



In the so called PT— symmetric quantum mechanics [2] one introduces an 
additional requirement which may be mathematically interpreted as the as- 
sumption of the Hermiticity of the Hamiltonian in a suitable Krein space 
[niEZ]- Formally speaking, one just preselects an indefinite Krein-space met- 
ric V with the property = I (called, conventionally, "parity"). One then 
postulates the following indefinite-metric parallel of Eq. ( ITO|) . 



[H]^V = VH. 



(17) 



The main benefit of such an auxiliary assumption is seen in the possibility 
of the construction of a special and unique metric ©("^^T") = CV called CVT 
metric. The new operator C with property = / is being interpreted as a 
charge [28] . 

Our present toy-model Hamiltonian matrices H = H^^\a, z) satisfy 
Eq. (11 7p with the Hamiltonian-independent matrices of parity V = 'p(^) con- 
taining just the unit elements along the secondary diagonal (i.e., Vm,n = 1 iff 
m + n = + 1 while Vmn = otherwise). Thus, at A^ = 2 one could treat 

as an indefinite limit of 6 [H^'^\a)] with fc(^) while 
and k'^^^nS^'^ — )■ 1. From the factorization requirement 0('^'^'^) = CV we may 
eliminate the complex charge 



U V 

y z 



18) 



The A^ 



and get v = y = k and z = u* = fce*"^ cos /3 (cf. Eq. 
version of condition = I requires not only that 7 = 'j^'^'^'^^ = 7r/2 [i.e., 
that a = cos/3) but also that /3 = is such that sin;^^^^^) = l/k. The 

resulting metric is unique and we have the unique charge 



c(^^^) = k ■ 



— ia 1 
1 ia 



(19) 



1 T 



The comparison of this formula with the general specification of an observable 
(fT6|) reveals that the charge is an observable in which D = b = c = and 
g = — cos/3/ sin/? = — a/A;2 — i. 



(k,m,r,h) 



6 The first nontrivial = 4 model 

6.1 The complete set of pseudometrics 

The general four-parametric Hermitian = 4 candidates for the metric (= 
"pseudometrics") may be obtained, most easily, from Eq. ffTOl) again. These 
solutions appear symmetric with respect to their second diagonal. They may 
be written in the closed four-parametric form of matrix G ^H^^^ (a, z) 

k m — ikw W* Z* 

m + ikw r h — i {kw + ra) W* 

W h + i {kw + ra) r m — ikw 

Z W 171 + ikw k 

where we introduced function w = w{z, a) = 3^a (this quantity must be 
positive and, for z > 0, larger than a) and where we abbreviated 

W = W{k, m, r) = —w^k + r — k — kwa + i {wm + ma) , 

Z = Z{k, m, r, h) = ma^ 



(20) 



w'^m 



m + h — i [kw — ka — kwa^ — rw + w^k^ . 

This result seems to indicate that the above-noticed full-matrix structure of 
the N = 2 metrics will survive the transition to any dimension. Our present 
choice of the discrete but strictly /oca/ non-Hermitian interactions V{xj) only 
seems to admit the non-band-matrix, strongly nonlocal forms of all of the 
Hermitizing metrics. 



6.2 The condition of positivity 

Although an immediate correspondence between the above-displayed general 
= 4 pseudometric with the most common and usual metric = / seems 
out of question, we may still select, for illustration purposes, the unit main 
diagonal in formula (120]) . k = r = 1. Once we also put m = h = (lead- 
ing to the simplified elements W{1, 0, 1) = —w {w + a) and Z{1, 0, 1, 0) = 
i{a + wa^ — w^)), we may immediately (and non- numerically) test the posi- 
tivity of the matrix. 

This test may proceed via an immediate evaluation of the four eigenvalues 

a, z)'j g ^ in the respective closed forms 



9^ of the candidate matrix B ^H^'^^ ( 



and 



= 1 + ^ {w-a^w + w^)±l ^/A+ 



= 1 - ^ {w - a^w + w"^) >/A- 



w 




Figure 15: The thick- line upper boundaries a = a{a) and (3 = P{a) of the 
respective positivity domains of the respective eigenvalues 9^ and 9g of the 
eligible metric 6 [H^^\a, z)~\ q ^ 

with the common discriminant 

= m;*^ + (2 - 2 a^) + (±8 + 4 a) + (5 ± 8 a + 6 + a^) 

+ {4a + Aa^)w + Aa^ . 

These formulae enable us to specify the parametric domain of the necessary 
positivity of the metric 9 [H^'^\a, z)~\ ^^q^q-^ by elementary means. 

The results are sampled in Fig. [T5]where we displayed the part of the {z, a) 
plane where 9^ > so that just the two eigenvalues may get non-positive 
there. The inspection of this picture enables us to conclude that these eigen- 
values remain positive (so that the general pseudometric 6 [H^^\a, z)~\ q ^ q-, 
becomes tractable as the positive definite metric) in the triangular domain 
n|Jf^' (assuming that z > 0) ot Q (assuming that z > 1). Such a simpli- 
fication follows from the fact that at the two sample values of 2; = and 
of ^ = 1, the auxiliary wiggly lines just "translate" the a— dependence into 
w— dependence of the change-of-sign boundary since the auxiliary functions 
w = w{z, a) = 3^a (which enter the above closed formulae as abbreviations) 
are linear in a. 

7 The N = 6 model 

For the not too large matrix dimensions = 2K the technique of the con- 
struction of the general pseudometrics via the solution of the linear al- 
gebraic Eq. f lTOj) remains feasible and straightforward. With the growth of 
N the only difficulty emerges in connection with the printed presentation of 
the resulting multiparametric set of pseudometrics. For this reason it makes 
sense to find a sufficiently representative set of a few key parameters. For 
each such choice, moreover, it becomes necessary to determine a boundary 
of the domain V of these key parameters, inside which the pseudometric in 

1 n 



question remains positive definite and, hence, eligible as a metric in a certain 
"optimal" physical Hilbert space 

In our present setting the task is slightly simplified by the fact that the 
most relevant values of our first, "kinematical input" parameter a may be 
expected small and admitting, therefore, the use of perturbation approxi- 
mations. Secondly, we may take the variability of the second, "dynamical 
input" parameter z for granted. Nevertheless, as long as the choice of this 
exponent will be mostly dictated by applications (controlhng, e.g., the loss of 
stability of certain most fragile levels), its choice proves not too essential in 
our present methodical considerations. We shall often consider the "discrete 
square-well value" [29] z = as a sufficiently instructive and generic option. 

7.1 The metric at z = 

Even when we set our "methodically redundant" exponent z equal to zero, 
we cannot parallel Eq. (120 p and display the whole six-parametric matrix G 
resulting from the computer-assisted non- numerical solution of Eq. ffTOl) . For 
this reason we further employed the simplification used in subsection 16.21 and 
demanded that all the elements of the main diagonal of our special G are 
chosen equal to one. Although the resulting pseudometric still did not fit 
in the printed page, the separate matrix elements do and remain sufficiently 
compact, 

G2,i = G6,5 = m + ia , G3,2 = G5,4 = 4 mc? + m + d + 2ia , 

©4,3 = d + A mc? + m + r + 3 ia , 

©3,1 = ©6,4 = ~2 o?+2 ima , G4,2 = Gs^s = —6 o?—i (—8 ma^ — 4 ma — 2 da) , 

©4,1 = ©6,3 = d-i{Aa^ - a) , ©5,2 = r + d-i (-2 a + 4 a^) , 

Gs^i = G6,2 = —4 a^ — i (—2 ma — 8 ma^ — 2 da) , Ge,! = r + ia . 

The complex conjugates of these elements form the upper-triangle part of 
G = G^^ of course. 
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Figure 16: The spectrum of the simplified metric (1211) . 
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In a continuation of the search of parallels with the results of subsection 



16.21 we may further set m 
sufficiently compact matrix 



d 



and obtain G = Qi^Ha) as the 
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—ia 


-2a2 


i (4a^ — a) 


-4a2 


— ia 


ia 


1 


—2 ia 


-6a2 


i (4a3 - 2a) 


-4a2 


-2a2 


2 ia 


1 


—3 ia 


-6a2 


i (4a^ — a 


i{a — 4a^) 


-6a2 


3 ia 


1 


—2 ia 


-2a2 


-A a' 


i (2a - 4a3) 


-6a2 


2 ia 
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— ia 


ia 


-4a2 


i (a - 4a3) 


-2a2 


ia 


1 



(21) 

The graph of its eigenvalues is displayed in Fig. [161 showing that such a 
matrix may serve as a metric in Hilbert space iff a G {— (3^^\ f3^^^) with 
^ 0.2718445. 




Figure 17: The spectrum of the perturbed metric ( l2Ti) with m = 1/10. 



The inspection of Fig. [16] reveals that in the domain of very small as 
the a— dependence of the eigenvalues 6j{a) of the metric 9o^''(a) > may 
be very well approximated by the linear functions, 6j (a) f» 1 + hj a. The 
determination of the coefficients hj is still easy since we may Taylor-expand 
the metric 



e(a) = J + 
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C (a^) (22) 



and arrive at the virtually trivial leading-order secular equation 



h^j -26h^ + 181 /i| - 225 = {h^j -2h]-llhj + 15) [h] + 2h]- 11 /ij - 15) = 



and roots = ^3. 102940862, /i2,5 = ±1.256913500, /i3,4 = ±3.846027361. 
The practical use of such a result may be twofold: 

• we may assure ourselves that the decrease of the lowest eigenvalue 
^jmini^) ~ 1 ^ 3.846 a guarantees that the maximal admissible value of 
a is /3 ^ 1/3.846027361 ^ 0.2600 plus/minus 0{a'^) corrections; this 
estimate seems fully consistent with the above-derived exact value of 
I3^^'> ^ 0.2718445; 

• alternatively, we may restrict ourselves to a much smaller subinterval of 
a in which the exact minimal eigenvalue Oj^.^ (a) remains safely positive, 
say, on the grounds of a variational estimate of corrections (keeping also 
in mind that maxa^ ^ 0.06). 

The formalism of perturbation expansions may be recalled and used to reflect 
the influence of further parameters. The importance of the possible changes 
of the spectrum of the metric is sampled in Fig. [17] where the differences 
from Fig. [16] are all caused by the choice of m = 1/10. One witnesses the 
complete removal of the full zero-order degeneracy of the "unperturbed" 
6j = 1 + 0{a). Just a partial removal of the degeneracy may result from the 
alternative choices, e.g., of r 7^ which just splits the six-times degenerate 
a = eigenvalue 9j in two three-times degenerate descendants. 

8 Extrapolations and models with > 8 

At any dimension N = 2K the dynamical contents of our toy models is con- 
trolled by the Hamiltonian (which varies with the "kinematical" parameter 
a e {—a^^\a^^^), "dynamical" parameter z e (— cxd, cxd) and real spectral 
locus {ej{a, z)}) andhj the Hermitian, positive-definite metric 9*^^^ specified 
by an A^— plet of parameters {k, m, . . .} G P. 

Naturally, practical implementations of such a recipe will require also a 
determination of the whole positivity domain V, or of its suitable non-empty 
subdomain at least. This task will certainly be facilitated by the smallness 
of a in practice. In our constructions performed at iV < 6 we saw, moreover, 
that many of the closed small— a formulae might be potentially extrapolated 
to an arbitrary dimension = 2K. 

8.1 The iV > 8 metrics 

One of the most promising keys to the extension of at least some of our above- 
listed low- dimensional results to all A^ has been found in the compact form of 
the matrix elements of the pseudometrics at A^ < 6. Indeed, we succeeded in 
transmuting these small— elements into larger— ansatzs and found out 
that such a method of construction appeared very efficient and quick. 

One of the main shortcomings of such a recipe lies in the enormous growth 
of the size of the formulae caused, mainly, by the linear growth of the number 
of free parameters with the increasing dimension A^ = 2K. This means 
that in the computer-assisted environment we may still deduce the form 



oo 



of the A^— parametric pseudometrics (from Eq. (fTO!) ) but the resulting 
extrapolations of the above-displayed compact K = 1, K = 2 and K = 3 
formulae cease to be compact. Thus, for presentation purposes, the majority 
of the available N > 8 results still has to be compactified and transformed, 
typically, into a graph or numerical table. 

During the computer-assisted N > 8 constructions themselves, the most 
difficult obstacle has been found, as we already indicated, in the necessary 
specification of the boundaries of a non-empty metric-positivity (sub) domain 
V^^K Fortunately, the very natural assumption of smallness of a almost 
trivialized the problem at < 6 and proved also helpful at the higher 
dimensions. The point was that the advantage of the practical negligibility 
of shortened the solution of Eq. (fTOj) . implying the - unexpected - easiness 
of the extrapolation of the matrix elements of pseudometrics to all A^. 

As we have already noted, the complete, multiparametric matrices G can 
hardly be displayed, in print, even at K = 3. For this reason, let us restrict 
our attention to certain special subsets of metrics and present the theory 
"via examples". Firstly, let us skip the questions of energies (discussed, at 
length, in the preceding sections after all) and study, from now on, just the 
simplest possible toy-model dynamics with z = 0. 

Secondly, let us circumvent the (difficult) problem of the determination 
of the exact "exceptional point" boundaries dV^^^ (redirecting the interested 
readers, say, to our dedicated study [30]) and let us select just a single repre- 
sentative point (i.e., multiindex /z := (fc,m, ...)—)■ {kQ,mQ, . . .) := /ig) which 
lies safely inside the hypothetical parametric domain T)^'^\ 

Thirdly, the guarantee of the latter requirement /iq € T)^^'> will be made 
easy via the extension of our previous A^ < 6 experience to all A^ and by the 
selection of the main-diagonal elements of the candidate for the metric equal 
strictly to one, Qn,n = I, n = 1,2, . . . , N (i.e., ko = 1, . . .), with all of the 
other parameters staying "sufficiently small" (i.e., |mo| -C 1, . . .). 

8.2 The > 8 metrics at z = and at small a 

In a continuation of our simplified presentation of the extrapolation ideas let 
us now select the dimension A^ = 8 and set all of the "small" parameters 
in a general pseudometric 9 strictly equal to zero (i.e., rriQ = 0, . . .). Imme- 
diately, the explicit solution of Eq. ffTOl) will generate the pseudometric Oq^^ 
with matrix elements which become predictable, by extrapolation, from their 
above-displayed A^ = 6 predecessors. 

On this basis we verified, numerically, that also the a— dependence of the 
spectrum {9j} of the new pseudometric (a) remains very similar to the 
one sampled in Fig. [16] above. We deduced, extrapolated and also re- verified 
at A^ = 10 that at any A^ = 2K, one may expect and conjecture to see the 
K—p\et of curves 6j{a) which are moving quickly up with the growth of \a\, 
complemented by another 7^— plet of curves Oj{a) which are moving quickly 
down with the growth of \a\. 

At this moment one can recollect the A^ = 6 discussion of Fig. [16] and 
formulate the following two questions concerning the determination of the 
pseudometric-positivity interval of a G {— (3^'^\ (3^^^) at the general dimen- 



sion: 



• do the left and right intersections ±(3^^^ G dV^^^ of the lowest eigen- 
value curve Oj^^^{a) with the a— axis leave a non-empty and/or suffi- 
ciently large space for the variability of a at large = 27^^? 

• would the linear approximation Oj^^.^^{a) = 1 + hj^.^a of the min- 
imal eigenvalues provide a sufficiently reliable estimate of the exact 
exception-point values of l3^^^7 

Whenever both of the answers happen to be positive, we may proceed 
further, accept the above-introduced perturbation-approximation philoso- 
phy and linearize our pseudometric. Thus, we returned to Eq. f fTOj) . con- 

(s) 

structed the linearized matrix 0q and obtained it in the following, very 
regular sparse-matrix form 
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We see that the extrapolation of this leading-order form of the pseudomet- 
ric to any = 2K is truly trivial. Also the construction of its linear- 
approximation eigenvalues Oj{a) = l+hj a remains feasible at any dimension. 
For illustration let us just select N = 8 and display the related leading-order 
secular equation for the coefficients, 

h^j - 70 h^j + 1487 /i^ - 9139 /i| + 11025 = . 

It may again be factorized yielding the two parallel rules 

h^. T2h^j- 33 h^j ± 47 hj + 105 = 

and roots 1.259204635, 2.752948888, 5.256297172 and 5.762552919 of both 
signs which determine the leading-order eight-line-crossing small— a part of 
the N = 8 analogue of the A^ = 6 spectra of Fig. [161 

We may add that the linear-extrapolation prediction 1/5.76255 ~ 0.1735 
of the exceptional-point value as obtained from the maximal root again com- 
pares very well with the exact numerical value of /3^^^ ~ 0.1683983. Thus, the 
linear approximation leads to the error ~ 0.005 which happens to be much 
smaller than the rough estimate ~ [f3^^^~\ ~ 0.03 of the typical magnitude 
of the second order correction. 



0/( 



8.3 Discussion: the > 8 discrete approximants and 
an efficient spectral design 

Our present proposal of working with discrete, matrix quantum models of 
finite dimension has been initiated by certain formal difficulties encountered, 
in the literature, during the attempted assignments of a nontrivial Hilbert- 
space metric to a given, non-Hermitian differential Hamiltonian operator 
possessing the strictly real spectrum. Typically, these difficulties are being 
circumvented via additional assumptions (cf. our comments above, especially 
in section [2]). 

During the very initial step of our analysis (viz., during the replacement of 
Eq. dl]) by Eq. using a finite number of grid points N = 2K) an elementary 
parameter a emerged and specified a fixed length in our models. In compar- 
ison with the differential-operator predecessors of our models, a new quality 
emerged since the spectra of energies became "fragile", i.e., complex beyond 
a certain critical kinematics-related size-parameter acHUcai = cn^^^ > 0. 

From the point of view of the variability of dynamics our initial choice 
of the one-parametric discrete versions of the very special potential of the 
imaginary cubic oscillator (as made in section [3]) did not prove too satisfac- 
tory. Fortunately, we managed to shift the role of the most unstable state 
to optional excitations by means of a re-scaling of the potential based on an 
introduction of another, "dynamical" exponent-parameter z G M. 

Undoubtedly, the latter trick made the structure (and, in particular, the 
"topological menu" ) of the energy levels "universal" in the sense illustrated, 
at = 8, in Table [2j At the same time, the mechanism of the changes 
of the topological structure of the spectral loci ej{a) remained transparent 
and tractable, schematically, as an up and down "motion" of the two (de- 
formed) circles. In this context the information compressed in Table [2] was 
complemented, graphically, by a series of illustrative pictures. All of these 
observations will find their strict analogues at any integer N = 2K. 

Another merit of our toy-model simulations of dynamics may be seen in 
the related feasibility of the constructions and in the extrapolation-friendliness 
of the multiparametric matrices of the pseudometrics. This implies, certainly, 
the rarely encountered and equally rarely employed freedom of the control 
of dynamics by the metric and of the related, rarely emphasized p3] theo- 
retical possibility of the ad hoc modifications of the classes of the additional 
observables Ai, A2, . . .. 

An important further merit of our present models has been found in the 
feasibility of the computer-assisted generation and extrapolations of impor- 
tant formulae (as well as of their graphical pendants and topological inter- 
pretations) to all N. As an immediate consequence one must appreciate, 
among others, the facilitated estimates of the ranges of the positivity of the 
pseudometrics, or the facilitation of the applicability of the linear-algebraic 
and perturbation-expansion techniques. 

Many of these ideas may find further applications. At the same time, we 
would expect that in the spirit, say, of Refs. [311 132], the next-step develop- 
ments of the subject should be aimed at the simulations of dynamics which 
replace the local potentials V{x) by some slightly non-local generalizations. 

OK 



Also in this respect, the present technique of discretization may be expected 
to lower many purely technical obstacles. 

9 Summary 

Our first tests of the idea of discretization proved disappointing. We observed 
that the growing-dimension series ([7]) of the simple-minded one-parametric 
descendants of the popular differential imaginary cubic oscillator do not of- 
fer a sufficiently rich variability of the coupling-dependence of the energy 
spectra. Fortunately, the merely slightly more sophisticated and re-scaled 
two-parametric choice ([8]) of the model has been shown to offer a flexibility 
of the spectral loci which covers a broad menu of alternative mechanisms of 
the phenomenologically interesting tunability. 

In our two-parametric model the breakdown of stability of the system was 
shown to be caused by the spontaneous, dynamically controlled complexifi- 
cation which could be designed as destroying the reality and stability of any 
pre-selected pair of neighboring bound states. 

The price to be paid for such a highly welcome universality of the model 
lies in the necessity of a rather complicated construction of the physical ad hoc 
metrics G. In this setting we took the advantage of the efficiency of the direct 
solution of the Dieudonne's Eq. ffTOj) as reported in Ref. |31]. We discussed the 
related idea of reduction of the ambiguity of the metric via the requirement 
of its maximal computational friendliness, i.e., of its sparse-matrix structure. 
We tested this possibility and arrived at encouraging results here. 

One of the main formal advantages of our present class of models may 
be seen in the possibility of its detailed study at the smallest dimensions 
followed by the formulation of ansatzs an by their tests and successful trial- 
and-error extrapolations to arbitrary dimensions. A particularly efficient 
application of such a strategy has been found in the context of perturbation- 
series constructions where we made use of the fact that due to its Runge- 
Kutta-approximation origin, the parameter a may be truly considered very 
small in practice. 

Our analysis also confirmed expectations that for the real exponents z, 
the metric operators can never be diagonal or banded matrices. This seems 
to be a characteristic consequence of the choice of a local form of the in- 
teraction for which the inner products remain "long-ranged" in the sense 
explained, in the context of scattering theory, by Jones [20] . Thus, in accord 
with our commentary in [19] we believe that the requirement of the unitar- 
ity of the scattering implies the necessity of introduction of at least small 
non-localities in the interaction. In opposite direction we expect that our 
present complex local-like Hamiltonians will preserve the analogy with their 
differential-operator-like N ^ 1 limits. In particular, we are persuaded (and 
would like to conjecture) that these models of dynamics will never admit the 
existence of a unitary and causal version of the scattering, not even in the 
"short-range" dynamical regime with very negative exponents z <^0. 
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